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1 Introduction 



It is well known that different types of topological objects may have been formed in the early 
universe after Planck time by the vacuum phase transition 012]. Depending on the topology 
of the vacuum manifold these are domain walls, strings, monopoles and textures. Among them, 
cosmic strings and monopoles seem to be the best candidate to be observed. A global monopolc 
is a spherical heavy object formed in the phase transition of a system composed by a self-coupling 
Goldstone field, whose original global symmetry is spontaneously broken. The matter fields play 
the role of an order parameter which outside the monopole's core acquires a non- vanishing value. 
The global monopole was first introduced by Sokolov and Starobinsky [3^. A few years later, the 
gravitational effects of of the global monopole were considered in Ref. 0], where a solution is 
presented which describes a global monopole at large radial distances. The gravitational effects 
produced by this object may be approximated by a solid angle deficit in the (3+l)-dimensional 
spacetime. 

The nontrivial properties of the vacuum are among the most important predictions of quan- 
tum field theory. These properties are manifested in the response of the vacuum to the external 
electromagnetic and gravitational fields. In particular, the explicit calculations of the vacuum 
polarization caused by particular external fields have played an important role in the devel- 
opment of quantum field theory. The quantum effects due to the point-like global monopole 
spacetime on the matter fields have been considered for massless scalar [3] and fermionic [5] 
fields, respectively. In order to develop this analysis, the scalar respectively spinor Green func- 
tions in this background have been obtained. The influence of the non-zero temperature on these 
polarization effects has been considered in [Jj for scalar and fermionic fields. Moreover, the cal- 
culation of quantum effects on massless scalar field in a higher dimensional global monopole 
spacetime has also been developed in [S]. The combined vacuum polarization effects by the non- 
trivial geometry of a global monopole and boundary conditions imposed on the matter fields are 
investigated as well. In this direction, the total Casimir energy associated with massive scalar 
field inside a spherical region in the global monopole background have been analyzed in Refs. 
|3 E3 by using the zeta function regularization procedure. Scalar Casimir densities induced 
by spherical boundaries have been calculated in I12j to higher dimensional global monopole 
spacetime by making use of the generalized Abel-Plana summation formula |131 HI] . More re- 
cently, using also this formalism, a similar analysis for spinor fields with MIT bag boundary 
conditions has been developed in ^3 ^] . 

Many of treatments of quantum fields around a global monopole deal mainly with the case 
of the idealized point-like monopole geometry. However, the realistic global monopole has a 
characteristic core radius determined by the symmetry braking scale at which the monopole is 
formed. A simplified model for the monopole core where the region inside the core is described 
by the de Sitter geometry is presented in |17j . The vacuum polarization effects due to a massless 
scalar field in the region outside the core of this model are investigated in Ref. JHj- in particular, 
it has been shown that long-range effects can take place due to the non-trivial core structure. In 
the present paper we will analyze the effects of global monopole core on properties of the quantum 
vacuum for the general spherically symmetric static model with a core of finite radius. The most 
important quantities characterizing these properties are the vacuum expectation values of the 
field square and the energy-momentum tensor. Though the corresponding operators are local, 
due to the global nature of the vacuum, the vacuum expectation values describe the global 
properties of the bulk and carry an important information about the structure of the defect 
core. In addition to describing the physical structure of the quantum field at a given point, the 
energy-momentum tensor acts as the source of gravity in the Einstein equations. It therefore 
plays an important role in modelling a self-consistent dynamics involving the gravitational field. 
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As the first step for the investigation of vacuum densities we evaluate the positive frequency 
Wightman function for a massive scalar field with general curvature coupling parameter. This 
function gives comprehensive insight into vacuum fluctuations and determines the response of a 
particle detector of the Unruh-DeWitt type moving in the global monopole bulk. The problem 
under consideration is also of separate interest as an example with gravitational and boundary- 
induced polarizations of the vacuum, where all calculations can be performed in a closed form. 
The corresponding results specify the conditions under which we can ignore the details of the 
interior structure and approximate the effect of the global monopole by the idealized model. 

The paper is organized as follows. In Section [2] we consider the Wightman function in 
the exterior of the global monopole for the general structure of the core assuming that the 
components of the metric tensor and their derivatives are continuous at the transition surface 
between the core and the exterior. By using this function, in Section 01 we investigate the 
vacuum expectation values of the field square and the energy-momentum tensor. The Section 
E] is devoted to the generalization of the corresponding results when an additional surface shell 
is present on the bounding surface between the core and the exterior. As an illustration of the 
general results, in Section [3 we consider the flower-pot model with the Minkowskian geometry 
inside the core. For this model the vacuum expectation values inside the core are investigated as 
well. In Section we present our concluding remarks. In Appendix we show that the formulae 
obtained in the paper for the core induced parts are also valid in the case when bound states 
are present. 

2 Wightman function 

We consider a model of (D + l)-dimensional global monopole with a core of radius a in which 
the spacetime is described by two distinct metric tensors in the regions outside and inside the 
core. In the hyper spherical polar coordinates (r, (f>) = (r, 0\, 02, ■ ■ ■ Q n , <^)> n = D — 2, the 
corresponding line element in the exterior region r > a has the form 

ds 2 = dt 2 - dr 2 - a 2 r 2 dfl 2 D , (1) 

where d£l 2 D is the line element on the surface of the unit sphere in D-dimensional Euclidean 
space, the parameter a is smaller than unity and is related to the symmetry breaking energy 
scale in the theory. The solid angle corresponding to Eq. (^Q) is a S D with S D = 2ir D l 2 /T{D/2) 
being the total area of the surface of the unit sphere in D-dimensional Euclidean space. This 
leads to the solid angle deficit (1 — a 2 ) So m the spacetime given by line element (0). It is 
of interest to note that the effective metric produced in superfluid 3 He — A by a monopole is 
described by the three dimensional version of the line element Q with the negative angle deficit, 
a > 1, which corresponds to the negative mass of the topological object [T^]. The quasiparticles 
in this model are chiral and massless fermions. We will assume that inside the core (region r < a) 
the spacetime geometry is regular and is described by the general static spherically symmetric 
line element 

ds 2 = e 2u ^dt 2 - e 2 < r Ur 2 - e 2w ^dQ 2 D , (2) 
where the functions u(r), v(r), w(r) are continuous at the core boundary: 

u(a) = v(a) = 0, w(a) = ln(aa). (3) 

Here we assume that there is no surface energy-momentum tensor located at r = a and, hence, 
the derivatives of these functions are continuous as well. The generalization to the case with 
an infinitely thin spherical shell at the boundary of two metrics will be discussed in section 
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Note that by introducing the new radial coordinate f = e w ^ with the core center at r = 0, the 
angular part of the line element (|2J) is written in the standard Minkowskian form. With this 
coordinate, in general, we will obtain non-standard angular part in the exterior line element (|T|). 
For the metric corresponding to line element (J2J) the nonzero components of the Ricci tensor are 
given by expressions (no summation over i, we adopt the convention of Birrell and Davies [2U] 
for the curvature tensor) 



vv + (n + 1) {w" + w' 2 — w'v')] , (4) 
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where the prime means the derivative with respect to the radial coordinate r and the indices 
i = 2, 3, . . . , D correspond to the coordinates 6±, 62, ■ ■ ■ , 4> respectively. The corresponding Ricci 
scalar has the form 

R = -2e~ 2v [v" + v' 2 - v'v' + n(n + l)w' 2 /2 

+(n + 1) (w" + w' 2 + w'v - w'v')] + n(n + l)e~ 2w . (5) 

Note that from the regularity of the interior geometry at the core center one has the conditions 
v(r),v(r) — > 0, and w(r) ~ lnf for f — * 0. In the region outside the core, r > a, for the nonzero 
components we have the standard expressions (no summation over i): 

1 2 1 <j 2 

R i = n -^r> R = < n + l )^pri ( 6 ) 

where i = 2, 3, . . . , D. For n = the spacetime outside the core is flat and coincides with D = 2 
cosmic string geometry. The influence of the non-trivial core structure for the cosmic string on 
a quantum scalar field has been considered in Refs. |21M22l 123] . In the discussion below we will 
assume that n > 0. 

In this paper we are interested in the vacuum polarization effects for a scalar field with general 
curvature coupling parameter £ propagating in the bulk described above. The corresponding 
field equation has the form 

(V.V* + m 2 + £R) tp = 0, (7) 

where Vj is the covariant derivative operator associated with line element (0) outside the core 
and with line element (J2J) inside the core. The values of the curvature coupling parameter £ = 0, 
and £ = £d with = (D — 1)/4D correspond to the most important special cases of minimally 
and conformally coupled scalar fields, respectively. As a first stage for the evaluation of the 
vacuum expectation values (VEVs) for the field square and the energy-momentum tensor we 
consider the positive frequency Wightman function (0\tp(x)(p(x')\0) , where |0) is the amplitude 
for the corresponding vacuum state. This function also determines the response of the Unruh- 
DeWitt type particle detector at a given state of motion (see, for instance, [20]). By expanding 
the field operator over eigenfunctions and using the commutation relations one can see that 

(0| 9 9(xMx')|0) = ^ V 9 Q (x)^(x / ), (8) 

a 

with {ip a (x), ip* a (x')} being a complete orthonormalized set of positive and negative frequency 
solutions to the field equation. The collective index a can contain both discrete and continuous 
components. In Eq. © it is assumed summation over discrete indices and integration over 
continuous indices. 
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Due to the symmetry of the problem under consideration the eigenfunctions can be presented 
in the form 



fa(x) = fi{r)Y(m k ;d,< 



-iujt 



1 = 0,1,2, 



where m& = (mo = I, mi, . . . , m n ), and m\,m2, ■ ■ ■ , m n are integers such that 
< m n _i < m n _ 2 <■■■ <mx<l, -m n -i <m n < m n -\, 



(9) 



(10) 



and Y (nik', $ , 4>) is the hyperspherical harmonic of degree I |24| . The equation for the radial 
function is obtained from the field equation Q and has the form 



//'(r) + [«' — v' + (D — l)w') fl(r) + e 2v [e 



m 



ZR-l(l + n)e- 2w ]Mr 



0. 
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In the region r > a described by the line element (^Q), the linearly independent solutions to this 
equation are r~ n / 2 J Ul (\r) and r~ n l 2 Y Vl (Ar) with A = \/u) 2 — m 2 , where J Vl { x ) an d ^ Ul (x) are 
the Bessel and Neumann functions with the order 



1 

n = - 
a 



n 



1 + +(l-a 2 )n(n + l)(e-^-i) 



(12) 



In the following consideration we will assume that v 2 is non-negative. This corresponds to the 
restriction on the values of the curvature coupling parameter for n > 0, given by the condition 

■2 



no 



4(n + 1) 



(13) 



This condition is satisfied by the minimally coupled field for all values a and by the conformally 
coupled field for a ^ D — 1. The solution of the radial equation ((TTj) in the region r < a regular 
at the origin we will denote by Ri(r, A). From Eq. (|11|) it follows that near the core center this 
solution behaves as r l . Note that the parameter A enters in the radial equation in the form A 2 . 
As a result the regular solution can be chosen in such a way that Ri{r, — A) = const • Ri(r, A). 
Now for the radial part of the eigenfunctions one has 



fl(r) 



Rl(r,\) 



for r < a 



-™/ 2 [AtJ Ul (Ar) + B{Y Vl (Ar)] for r > 



(14) 



where the coefficients A[ and B\ are determined by the conditions of continuity of the radial 
function and its derivative at r = a. From these conditions we find the following expressions for 
these coefficients 



A l = |a"/ 2 ^(a,A)^(Aa), B, = —a^R^a, \)J Ul (\a). 
Here and in what follows for a cylinder function F(z) we use the notation 



F{z) = zF'(z) 



n R[{a,z/a) 
2 +a R l (a,z/a) 



F(z), 



(15) 



(16) 



where i?^(a,A) = dRi(r, X)/dr\ r=a . Note that due to our choice of the function Ri(r, A), the 
logarithmic derivative in formula (|16|) is an even function on z. Hence, in the region r > a the 
radial part of the eigenfunctions has the form 



Mr) 



n/2 



2 r n/2 



Ri(a, X)g Ut (Xa, Ar) 



(17) 
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where the notation 

g Ul (Xa,Xr) = J U[ (Xr)Y U[ (Xa) - J U[ (Xa)Y Ul (Xr). (18) 

is introduced. 

For the eigenfunctions we have the following orthonormalization condition 

dV^99 00 Va(x)<p* af (x) = 6 -^-, (19) 

where 6 aa > is understood as the Kronecker symbol for discrete indices and as the Dirac delta 
function for continuous ones. Substituting eigenfunctions Q, and using the relation 

cM |y(m fc ;tf,0)| 2 = 7V(m fc ) (20) 

(the explicit form for N{mf t ) is given in [21] and will not be necessary for the following consider- 
ation in this paper), the normalization condition is written in terms of the radial eigenfunctions 

^ir^MnmM^^ly (2D 

where r$ is the value of the radial coordinate r corresponding to the origin and g r is the radial 
part of the determinant g. Note that in general ro ^ (see, for instance, the special case of 
the flower-pot model in section [SJ. As the integral on the left is divergent for A' = A, the main 
contribution in the coincidence limit comes from large values r. By using the expression (|17|) 
for the radial part in the region r > a and replacing the Bessel and Neumann functions by the 
leading terms of their asymptotic expansions for large values of the argument, it can be seen 
that from (|21j) the following result is obtained: 

2rr 1 ~ D X 

° njR ' (a ' A) = 2 N< U72M UV2a \Y (22) 

■K 2 u)N(m k ) [J4(Aa) + Yj ( (\a)\ 

Having the normalized eigenfunctions, now we turn to the evaluation of the Wightman 
function by using the mode sum formula jSJ). Substituting eigenfunctions (|17|) and using the 
addition formula for the hyperspherical harmonics |24j 

\ - y(m fc ; , , 2l + n w / 2 , . . , 

\ N(m k ) Y{mk ^ A) = ^ Cl (C ° S0) ' (23) 

for the Wightman function in the region outside the monopole's core one obtains 

l—D °° 9/ , 

(oMxMx'M = ^Et^^ /2 (cos^) 

1=0 v ' 

XdX g Vl (Aa, Xr)g Ul (Xa, Xr') iy /j7+^ (t >_ t) , , 

^WT^ J2(Aa) + ^(Aa) ' 1 > 

In formula S D = 2tt d ^/T(D/2) is the total area of the surface of the unit sphere in D- 
dimensional space, Cp(x) is the Gegenbauer or ultraspherical polynomial of degree p and order 
q, and 9 is the angle between directions ($,0) and ($', <j>'). Let us denote by (0 m \p(x)ip(x')\0 m ) 
the positive frequency Wightman function for the geometry of the idealized point-like global 
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monopole described by the line element (^Q) for all values of the radial coordinate. This function 
can be presented in the form jllj 

l—D 00 9/ 1 



00 \ W\ 2 +m 2 (t'-t) 

d\—===-J Vl (\r)J Vl {Xr'). (25) 
v X + m z 

In order to investigate the part induced by the non-trivial core structure, we consider the dif- 
ference 

(<p(x)<p(a/)) c = <0|¥>(xMz')|0) - (0 m \<p(x)<p(x')\0 m ). (26) 



Using formulae (|24fl . (|25jl and the relation 

9u (Xa Xr) g (_Xa Xr>) _ l£ J^±H^(Xr)H^(Xr'), (27) 

J2(Aa) + v(Aa) 2 ^ Hif'(Xa) 

with Hy S \x), s = 1, 2 being the Hankel functions, the core induced part in the Wightman 
function is presented in the form 



Z=0 



e 

x - 



7, lAaJ Hff(Xr)H$(Xr'). (28) 
v^T^ ^W(Aa) ' ' 



Now we rotate the integration contour in the complex plane A by the angle tt/2 for s = 1 and by 
the angle —tt/2 for s = 2. By using the property that the logarithmic derivative of the function 
i?/(r, A) in formula is an even function on z, we can see that the integrals over the segments 
(0, im) and (0, —im) of the imaginary axis cancel out. As a result, after introducing the modified 
Bessel functions, the core induced part can be presented in the form 



pi—Li ^ 21 -\- n / f°° 
(<p(x)(p(x')) c = 5— ^ 7 , C" /2 (cos6>) / dzz-^ 



K Ul (zr)K Ul {zr') 
V ' z 2 — m 2 



cosh 



Vz 2 -m 2 (t' -t) . (29) 



z 2 — m 2 

Here and below the tilted notation for the modified Bessel functions is defined as 

F(z) = zF'(z)-K l (a,z)F(z), (30) 

with 

, , n R',(a,ze m / 2 /a) . . 

The VEVs in the bulk of the idealized point-like global monopole are well-investigated in liter- 
ature (see, for instance, 0-^2] and references therein) and in the discussion below we will be 
mainly concerned with the part induced by the non-trivial core structure. As we see from (|29|) . 
all information about the inner structure of the global monopole is contained in the logarith- 
mic derivative of the interior radial function in formula In deriving formula (|29|) we have 
assumed that there are no bound states for which A is purely imaginary. In appendix we show 
that this formula is also valid in the case when bound states are present. 
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3 Vacuum expectation values outside the monopole core 



The VEV of the field square is obtained by computing the Wightman function in the coincidence 
limit x' — ► x. In this limit expression (|24j) gives a divergent result and some renormalization 
procedure is needed. Outside the monopole core the local geometry is the same as that for a 
point-like global monopole. Hence, in the region r > a the renormalization procedure for the 
local characteristics of the vacuum, such as the field square and the energy-momentum tensor, is 
the same as for the point-like global monopole geometry. This procedure is discussed in a number 
of papers (see |S] - [E])- For the renormalization we must subtract the corresponding DeWitt- 
Schwinger expansion involving the terms up to order D. For a massless field the renormalized 
value of the field square has the structure (</? 2 )m,ren = [A + B ln(/xr)] /r D ~ x , where the coefficients 
A and B are functions on the parameters a and £ only and the arbitrary mass scale fi corresponds 
to the ambiguity in the renormalization procedure. For a spacetime of odd dimension B = 
and this ambiguity is absent. In general, it is not possible to obtain closed expression for the 
coefficients A and B. For small values 1 — a 1 approximate expressions are derived in Ref. [H] 
for D = 4 and D = 5. In this paper our main interest are the parts in the VEVs induced by the 
non-trivial core structure and below we will concentrate on these quantities. 

By using the formula for the Wightman function from the previous section, the VEV of the 
field square in the exterior region is presented in the form 



Taking into account the relation 



r n/2 m _ T(l + n) 

Cl U) - T(n)l\ ' {S2} 



for the part induced by the core we find 

, ■>,,, K n 

Here the factor 



(A = -^EA dzz Wa) PS . (33) 



D , = (2i + D -2) ^ D1)I > (34) 

is the degeneracy of each angular mode with given I. For a fixed I and large z the integrand 
contains the exponential factor e 2z ( a_r ) and the integral converges when r > a. For large values 
I, introducing a new integration variable y = z/v\ in the integral of Eq. Q33|) and using the 
uniform asymptotic expansions for the modified Bessel functions it can be seen that the 
both integral and sum are convergent for r > a and diverge at r = a. For the points near 
the sphere the part (|33|) behaves as l/(r — a) , where Pi is an integer which depends on the 
specific model of the core. For this parameter one has j3\ ^ D — 1. The exception is the case of 
the core model for which the leading term in the uniform asymptotic expansion of the function 
K Vl (za) for large values I vanishes. The latter takes place for the interior radial function with the 
asymptotic behavior lZi(a,lz/a) ~ —{l/a)y/l + z 2 for large I. For the case of a massless scalar 
the asymptotic behavior of the part (|33|) at large distances from the sphere can be obtained 
by introducing a new integration variable y = zr and expanding the integrand in terms of a/r. 
The leading contribution for the summand with a given / has an order (a/r) 2l/ ' +£>_1 [assuming 
that v\ 7^ and lZi(a,0) ^ ±i/j] and the main contribution comes from the I = term. Now 
comparing this with the part (92 2 ) m ,ren> we see that for uq > the VEV of the field square at 
large distances from the core is dominated by the part corresponding to the geometry of the 
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point-like global monopole. For the case v q = the ratio ((p 2 ) c / ( l P 2 )m,ren decays logarithmically 
and long-range effects of the monopole core appear similar to those for the geometry of a cosmic 
string \'22\ (see also the discussion in Ref. ^H] for the model with de Sitter spacetime inside 
the core). This case is realized by special values of the parameters satisfying the condition 
(1/cr 2 — 1)£ = — £d— i- For a massive field assuming that mr ^> 1, the main contribution into 
the integral over z in Eq. (|33|) comes from the lower limit and to the leading order one has 



' (mo) 



(35) 



with the exponentially suppressed VEV. 

Consider the limit uC 1 for a fixed value r. In accordance with Eq. © this corresponds to 
large values of the scalar curvature and, hence, to strong gravitational fields. To satisfy condition 
(|13j) we will assume that £ ^ 0. For £ > 0, from Eq. (|12j) one has m 3> 1, and after introducing 
in Eq. ()33|) a new integration variable y = z/vi, we can replace the modified Bessel function 
by their uniform asymptotic expansions for large values of the order. The main contribution 
to the sum over I comes from the summand with I = 0, and the core induced VEV (if 2 

suppressed by the factor exp 



is 



For £ = and a <C 1 for the terms 



-(2/a) v / n(n + l)£ln(r/a) 
with I 7^ one has v\ 3> 1 and the corresponding contribution is again exponentially small. For 
the summand with / = to the leading order over a we have v\ = n/2 and (ip 2 ) c ~ l/o D ~ l . 
Hence, we conclude that in the limit of strong gravitational fields the behavior of the VEV (</? 2 ) c 
is completely different for minimally and non-minimally coupled scalars. 

Now we turn to the investigation of the VEV of the energy-momentum tensor in the region 
r > a. Having the Wightman function and the VEV for the field square, these VEVs are 
evaluated on the base of the formula 



(0|r«|0) = lim di& k (Q\<p(x)ip{x')\0) + 



e - 1 ] a-kViV 



(0\<p 2 (x 



(36) 



Similar to the Wightman function, the components of the vacuum energy-momentum tensor can 
be presented in the decomposed form 



<0|T ifc |0) = <0 m |T iJfc |0 m > + (T ik ) 



(37) 



where (OmlTj^lOm) is the vacuum energy-momentum tensor for the geometry of a point-like global 
monopole and the part {Ti k ) c is induced by the core. In accordance with the problem symmetry 
both these tensors are diagonal. For massless fields the VEV of the energy-momentum tensor 
for the point-like global monopole geometry is investigated in Refs. |51 _ jH]- The corresponding 
renormalized components have the structure similar to that given for the field square: 



(Ti 



1 



ifc/m,rcn 



-D+l 



o (1) 



+ 'iik 



ln(//r) 



(38) 



where the coefficients q. 



ik ' 



a {2) 



depend only on the parameters a and £, and q, 



(2) 
ik 



for D 



being an even number. Substituting the expressions of the Wightman function and the VEV of 
the field square into formula (|36|) . for the part of the energy- momentum tensor induced by the 
non-trivial core structure one obtains (no summation over i) 



in 



27rr n S D ^ 
u 1=0 



D, 



dz z 



3 I Vl (za) F® [K Ul {zr)) 



K n (za) \fz 1 



r > a, 



(39) 
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where for a given function f(y) the notations 



[f(y)} 
f£ ] lf(y)} 

f$> lf(y)] 



(1-40 



i'Ii \ n jn\£H\, ( v l l + 4£-2(mr/y) > 2 

/ \y) - -f(y)f (y) + ) f 2 (y) 



f 2 (y) + lf(y)f'(y) - ( 1 + "' 2+ | n/2 1 f\y), 



(^-i)f' 2 (y)--f(y)f(y) + 
y 



y 



4f _ 1+ ^ 2 (l + l)+W2 



(40) 
(41) 



(n + 1)^ 



/ 2 (y), (42) 



are introduced with £ = 4(n + 1)£ — n and in Eq. (|42jl 
components (f3l?)) satisfy the continuity equation Vk{T^) c 
consideration takes the form 



= 2, 3, . . . , D. It can be seen that 
0, which for the geometry under 



4^ 



+ ( J D-l)((T 1 1 ) c -(r 2 2 ) c )=0. 



(43) 



The core induced part (T^) c are finite everywhere outside the core, r > a, and diverge on 
the core boundary. Near this boundary the main contribution comes from large values I and to 
find the corresponding asymptotic behavior we can use the uniform asymptotic expansions for 
the modified Bessel functions. To the leading order one finds (T t k ) c ~ l/(r- af 2 for the energy 
density and the azimuthal stress and 



D - 1 



r/a 



l)(7f) c , 



(44) 



with fa ^ D + 1. An exception is the special case of the core model for which the leading term 
in the uniform asymptotic expansion for the function K Ul (za) vanishes. For large distances from 
the core boundary, r 3> a, and for a massless scalar field the main contribution into the VEV 
(Tf) c comes from the / = summand. Under the assumptions vq ^ and H$(a,Q) ^ iz'o, the 
leading term of the corresponding asymptotic expansion behaves like (T^) c ~ (a/r) 2u ° +D+l . For 
a massive scalar field under the condition mr 3> 1, the main contribution into the integral over 
z in Eq. (|39|) comes from the lower limit and by using the asymptotic formulae for the function 



K Ul (zr) for large values of the argument, to the leading order one finds 



(Tl 



)c 



-(Ti) c 



(4£ - 1; 



0p m 3/2 e - 



2mr 



4r D-l/2g Da 



I n (ma) 



1=0 



K n (ma) 



(45) 



and the radial stress is suppressed by an additional factor 1/mr. 

Now let us consider the VEV of the energy-momentum tensor in the limit a<l for a fixed 
r > a. For £ > by the calculations similar to those given above for the field square, one finds 
that the core induced VEV are suppressed by the factor exp — (2/a)^Jn(n + 1)£ ln(r/c 

the vacuum stresses are strongly anisotropic: (T 1 1 ) C /(T 2 2 ) c ~ a. For a minimally coupled scalar, 
£ = 0, the leading term of the asymptotic expansion over a comes from the I = summand in 
Eq. ()39|) with v\ = n/2. This term behaves as a l ~ D . 



and 



4 Core with an infinitely thin shell 

The results considered in the previous section can be generalized to the models where an addi- 
tional infinitely thin spherical shell located at r = a is present with the surface energy-momentum 
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tensor rf. We denote by n l the normal to the shell normalized by the condition n; L n % = — 1, 
assuming that it points into the bulk on both sides. From the Israel matching conditions one 
has 

{K lk - Kh lk ] = torGr ik , (46) 

where the curly brackets denote summation over each side of the shell, h{ k = + n{n k is the 
induced metric on the shell, Ki k = ^/ttV r n s its extrinsic curvature and K = K\. For the region 
r ^ a one has rii = 5je v ^ and the non-zero components of the extrinsic curvature are given by 
the formulae 

K$ = -u'(r)e-< r \ K\ = -5{w' '{r) e - v{r \ i = 2, 3, . . . , r = a - . (47) 

The corresponding expressions for the region r a are obtained by taking u(r) = v(r) = 0, 
w(r) = In (or) and changing the signs for the components of the extrinsic curvature tensor. Now 
from the matching conditions (|46|) we find (no summation over i) 



u'(o-) = 87TG 

1 8ttG 



r i D ~ 2 J0 



2,3,..., (48) 



*"<»-) " a + —l T ° < 49 » 

where /'(a—) is understood in the sense lim r ^ a _o /'(?*)• ^ ne discontinuity of the functions u'(r) 
and w'(r) at r = a leads to the delta function term 

2 [u'(a-) + (D - 1) (w'(a-) - l/a)] 5{r - a) (50) 

in the Ricci scalar and, hence, in the equation for the radial eigenfunctions. Note that 
the expression in the square brackets is related to the surface energy-momentum tensor by the 
formula 

u'(o-) + (D-l) (w'(a-) - l/a) = i^r, (51) 

where r is the trace of the surface energy-momentum tensor. 

Due to the delta function term in the equation for the radial eigenfunctions, these functions 
have a discontinuity in their slope at r = a. The corresponding jump condition is obtained by 
integrating the equation 1)1 ljl through the point r = a: 

f[(a+)-fl(a-) = ^frf l (a). (52) 

Now the coefficients in the formulae (|14p for the eigenfunctions are determined by the continuity 
condition for the radial eigenfunctions and by the jump condition for their radial derivative. It 
can be seen that the corresponding eigenfunctions are given by the same formulae Q17|) and (|22j) 
with the new barred notation 



F(z) = zF'(z) 



n 16ttG£ R(a,X) 

- H -ar + a l ) ( 

2 D-l Ri(a,X) 



F(z). (53) 



Consequently the parts in the Wightman function, in the VEVs of the field square and the 
energy- momentum tensor induced by the core of the finite thickness, are given by formula (|29jl. 
(|3"3|) and where the tilted notation is defined by Eq. PU|) with the function 



n 16ttG£ _ ; i^(a,ze™/ 2 /a) 



Hi (a, z) = - + -j^-ar + a-f- . (54) 

2 D — 1 Ri(a, ze m ' z /a) 

The trace of the surface energy-momentum tensor in this expression is related to the components 
of the metric tensor inside the core by formula (|51|) . 
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5 Flower-pot model for global monopole 



As an application of the general results given above let us consider a simple example of the core 
model assuming that the spacetime inside the core is flat. The corresponding model for the 
cosmic string core was considered in Refs. |21M22ll2*3] and following these papers we will refer to 
this model as flower-pot model. Taking u(r) = v(r) = from the zero curvature condition one 
finds = r + const. The value of the constant here is found from the continuity condition 
for the function w(r) at the boundary which gives const = (a — l)a. Hence, the interior line 
element has the form 

ds 2 = dt 2 - dr 2 - [r + (a - l)a} 2 dn 2 D . (55) 

In terms of the radial coordinate r the origin is located at r = (1 — a)a. From the matching 
conditions ([48)1 . (|49jl we find the corresponding surface energy-momentum tensor with the non- 
zero components 

T -(±-l]^± r k -^T°5 k i-2 3 (56) 
T ° " [a l ) 8nGa' ^ ~ D — 1 1 ~ l 5b J 

The corresponding surface energy density is positive for the global monopole with a < 1. After 
this brief review, let us analyze for this model the influence of the monopole's core on the vacuum 
polarization effects. We will consider the exterior and interior regions separately. 



5.1 Exterior region 

In the region inside the core the radial eigenfunctions regular at the origin are the functions: 

^+n/2(Ar) 



Rl{r,X) = C l - 



rn/2 



(57) 



where f = r + (a — l)a is the standard Minkowskian radial coordinate, ^ f ^ aa. In 
appendix we show that in the flower-pot model no bound states exist. Note that for an interior 
Minkowskian observer the radius of the core is era. The normalization coefficient C\ is found 
from the condition (|22() : 



C 



2AJ, 



-2 

l+n/2 



(Xaa) 



■n 2 aoJi 



;N(m k ) [JS(Aa) + y2(Ao)]' 
with the barred notation for the cylindrical functions 



(58) 



F(z) = zF'(z) 



a a + z 



J 'l + n/2( Z(J ) 



and 



Or, 



Jl+n/2{zcr 

ra - 4£(ra + 1)] . 



F(z) 



(59) 



(60) 



Note that J U[ (\a) = for a = 1. Hence, the parts in the Wightman function, in the VEVs of 
the field square and the energy-momentum tensor due the non-trivial structure of the core in 
the flower-pot model, are given by formulae (|29"|) . and (|39"|) respectively, where the tilted 
notations for the modified Bessel functions are defined by ()3fl|) with the coefficient 



TZi(a, z) = a a + z 



r i+n/2( Za ) 



h+n/2{ 



za 



(61) 
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For a = 1 one has I Vl (z) = and as we could expect the VEVs vanish. Using the value for the 
standard integral involving the product of the functions K v given in Ref. [25], in the case of a 
massless scalar field the leading term for the asymptotic expansion over a/r can be presented 
in the form 

2 . v F(2u + 1/2)I> + 1/2) An fa\ 2u +D-l 



where 



A n 



na - 4£(ra + l)(a - 1) - 2av Q 



(62) 



(63) 



na - 4£(n + l)(a - 1) + 2av 

Note that for a minimally coupled scalar A n = and the presented leading term vanishes. In 
figure ^ we have plotted the dependence of the part in the VEV of the field square induced by 
the core as a function on the rescaled radial coordinate for minimally and conformally coupled 
D = 3 massless scalar fields in the flower-pot model with a = 0.5 

0.0008 



0.0006 



0.0004 



0.0002 



-0.0002 



1.5 2 2.5 3 3.5 4 



Figure 1: The expectation value a D ~ 1 (ip 2 ) c induced by the non-trivial core structure in the 
region outside the core for D = 3 massless scalar field as a function of r/a in flower-pot model 
with a = 0.5. The full/dashed curves correspond to minimally /conformally coupled scalars. 

Now let us analyze the VEV of the energy- momentum tensor given by Eq. (|39|) with the tilted 
notation given by (|30|) . (|61|) . For large distances from the core, r 3> a, the main contribution 
into the VEV of the energy-momentum tensor for a massless scalar field comes from the I = 
summand. Under the assumption ^ the leading terms of the asymptotic expansions have 
the form (no summation over i) 



2- 2l/ °a 1 - D A n 5f 
TtvoS D T 2 (vQ)a 



D+l 



r 



dzz 2 »° +2 F$[K„ (z)]. 



(64) 



The integrals in this formula can be evaluated using the value for the integrals involving the 
product of the functions K v given in Ref. [22j. As we see, for v§ > and for large distances 
from the sphere the vacuum energy-momentum tensor is dominated by the part corresponding 
to the point-like monopole. 

As it has been mentioned above on the core surface the VEVs diverge. For the region near 
the core the main contribution comes from large values of I. By using the uniform asymptotic 
expansions for the modified Bessel functions it can be seen that to the leading order (y? 2 ) c ~ 
(r — a) 2 ~ D and the components of the vacuum energy-momentum tensor behave as (r — a)~ D 
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for the energy density and the azimuthal stress and as (r — a) 1_D for the radial stress. Due to 
surface divergencies near the surface the total vacuum energy-momentum tensor is dominated 
by the parts induced by the finite thickness of the core. As an illustration, in figure |21 we have 
presented the dependence of the core-induced vacuum energy density as a function on the radial 
coordinate for D = 3 minimally and conformally coupled massless scalar fields in the flower-pot 
model with a = 0.5. 

0.00 4 1 , — , , , , , , , , , , , , , , , 



0.002 




-0.002 ■ 



1.2 



Figure 2: The expectation value of the energy density, a D+1 (Tq) c induced in the region outside 
the core for D = 3 massless scalar field as a function of r/a in flower-pot model with a = 0.5. 
The full/ dashed curves correspond to minimally /conformally coupled scalar s. 



5.2 Interior region 

Now let us consider the vacuum polarization effects inside the core for the flower-pot model. 
The corresponding eigenfunctions have the form given by Eq. © with fi(r) = Ri{r, A) and the 
function Ri(r, A) is defined by formula (|57j). Substituting the eigenfunctions into the mode sum 
formula for the corresponding Wightman function one finds 



TT 2 naSn ^ (fr') n l 2 

°° A\ AJ '+ n / 2(Aafl) Ji+n/2^r)J l+n/2 (Xr') iV ^? {t ,_ t) 

To find the renormalized VEVs of the field square and the energy-momentum tensor we need to 
evaluate the difference between this function and the corresponding function for the Minkowski 
bulk: 

{<p(x)<p(x')) suh = {0\ip(x)cp(x')\0) - {0 M \<p(x)<p(x')\0 M ). (66) 

The appropriate form for the Minkowskian part is obtained from Eq. (|25[) taking a = 1 and 
replacing r — ► f. By using the corresponding formula, for the subtracted Wightman function 
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one finds 



(tp(x)tp(a/)) 



sub 



21 + n 



ir 2 nS D (fr') n / 2 



E 



q" /2 (oos«) / 

J 



oo \ W\' 2 +m' 2 (t'-t) 



X^+n/2(Ar)^+n/2(Ar') 



j2(Aa) + y2( Aa ) 4 



(67) 



where the barred notation is defined by Eq. (|59|). The integral in this formula is slowly con- 
vergent and the integrand is highly oscillatory. In order to transform the expression for the 
subtracted Wightman function into more convenient form, we note that the following identity 
takes place 



1 C{H^ n/2 (az),J Ul (z)} 

2 ^ ' 

s=l,2 



1, 



where we have introduced the notation 

C {f(crz),g{z)} = zf(az)g'(z) - [a ff f{az) + zf'(az)] g(z). 
Note that in terms of this notation one has 

Ji+n/2(o-z)F(z) = C{J l+n/2 {az),F(z)}. 



(68) 



(69) 



(70) 



We add the left-hand side of Eq. (|68[) with z = Xa as a coefficient to the term 7r 2 /4 in the square 
brackets of Eq. (|67|) . After this replacement the term in the square brackets is written in the 
form 



J, 



l+n/2 



(az)/a 



? - E 



^ 2 2C{J l+n/2 (az),J Ul (z)} 
C{J l+n/2 (az),H$(z)} 



7T" 



C{H^ n/2 (<yz),J Vl (z)} 



(71) 



Note that both terms in the sum over s on the right of this relation are separately regular at 
the zeros of the function C{ J;_|_ n / 2 (cAa), J Ul (Xa)}. Substituting (|7T|) into formula (|oT|) we rotate 
the integration contour in the complex plane A by the angle tt/2 for s = 1 and by the angle 
— 7r/2 for s = 2. Under the condition r + f' + \t — 1'\ < 2aa the contribution from the semicircle 
with the radius tending to infinity vanishes. Note that as we consider the points inside the core 
this condition is satisfied in the coincidence limit. The integrals over the segments (0, im) and 
(0, —im) of the imaginary axis cancel out and after introducing the modified Bessel functions 
the subtracted Wightman function can be presented in the form 



{ip(x)ip(x')) suh 



_J_X^ 2l + n n/2, 

/ ~, ; — 7T7<-v; COS I 

TmS D ^ (ff'Y t2 

1=0 y ' 



00 , zUi(a,za) 
dz 



with the notation 



Ui(a,z) 



xll+n/2(zr)li +n / 2 {zr) cosh \J z 2 -m 2 (t' - t) 

1/a + C{I l+n/2 (az),K Ul (z)}C{K l+n/2 (az),I Ul (z)} 
C{Ii +n /2(vz),I»i (z)}C{Ii +n/2 (az),K Ul (z)} 



(72) 



(73) 
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For points away from the core boundary the integral is exponentially convergent in the coinci- 
dence limit and this formula is convenient for the calculation of the VEVs of the field square 
and the energy-momentum tensor. 

Having the subtracted Wightman function we can evaluate the renormalized VEV of the 
field square by taking the coincidence limit of the arguments in Eq. (|72[): 

^=-iK^ w * (74) 

So, by this result we can see that although the spacetime inside the core is Minkowski one, there 
exists a vacuum polarization induced by the non-trivial topology of the spacetime in the exterior 
region. In the limit r — > the main contribution into Eq. (|74|) comes from the I = summand 
with the leading term 

^ >™ ra - ^WrP) I 7F=^' (75) 

and at the monopole center the renormalized VEV of the field square acquires a non- vanishing 
regular contribution. For the points on the core surface the VEV given by Eq. (|74|1 diverges 
like \/{a — r) D ~ 2 . Now let us consider the limiting case a <C 1 under the fixed value aa which is 
the core radius for an internal Minkowskian observer. Introducing in Eq. (|74|) a new integration 
variable y = aaz and by making use of the uniform asymptotic expansions for the modified 
Bessel functions with the index v\, we can see that to the leading order 



{v >- " g D ' L dz tw^? iL^rr (76) 



where we have introduced the notation 



m + y - m + y^w/^iyn 

with 

? ? i = I (I + n) + n(n + 1)£, ? ?2 = 2^(n + 1) - n/2. (78) 

Hence, in the limit a — > for a fixed radius of the core, aa, the part in the renormalized VEV 
of the field square inside the core tends to the finite limiting value. For large values of the 
mass, assuming that m(aa — f) 3> 1, it can be seen that (<p 2 ) ven is suppressed by the factor 
e -2m(aa-r) _ j n f|g ure we have plotted the renormalized VEV (f 2 ) r en inside the core of the 
flower-pot model with a = 0.5 as a function of f/aa for minimally and conformally coupled 
massless scalars. Again we can observe that there exists a strong dependence of this quantity 
on the curvature coupling parameter. 

The renormalized VEV of the energy-momentum tensor is found by using the formula (|36j) 
with the subtracted Wightman functions. This leads to the following formula (no summation 
over i) 

where the functions F^ n , 2 [f(y)] are defined by relations (|4Uj ) -(|42j l with the replacement i>\ — ► 
I + n/2. At the core center the nonzero contribution to VEV (|79|) comes from the summands 
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Figure 3: The expectation value aP 1 ((/? 2 )ren inside the core for D = 3 massless scalar field as 



a function of r/aa in the flower-pot model with a 
minimally /confer mally coupled scalar s. 



0.5. The full/dashed curves correspond to 



with I = and I = 1 and one has 

(^o')ren : 



2 D n D / 2 + 1 T{D/2) J m 



00 z D+l dz 



4£ + 1 - 2- 



m 



U (a, za) + (4f - l)t/i(<r,zo) 



(Ti 



2 /ren 



2 D n D / 2 + 1 DT(D/2) J m ^ 



00 z D+1 dz 



£-2 l7o(<r, za) +£C/i (<r,za) 



(80) 



(81) 



Note that for the conformally coupled massless scalar at the center one has (TQ) ren = — Z?(T 1 1 ) reri . 
This can also be obtained directly from the zero trace condition. Near the core surface the com- 
ponents of the vacuum energy-momentum tensor behave as (a — r)~ D for the energy density and 
the azimuthal stress and as (a — r) l ~ D for the radial stress. As in the case of the field square, in 
the limit a — > for a fixed radius of the core radius act, the part in the vacuum energy- momentum 
tensor induced by the non-trivial core tends to the finite limiting value. This limiting value is 
obtained from formula (|T§)) by making the replacement Ui(cr,za) — ► Ui +n /2{zaa) / : I 2 +n ^ 2 (zaa). 

As in the case of the field square, for large values of the mass for the field quanta the VEV (T^) ien 
is exponentially suppressed by the factor e - 2m ( CTa - r ) . The dependence of the renormalized inte- 
rior vacuum energy density on the radial coordinate is presented in figure 0] for minimally and 
conformally coupled massless scalar field in D = 3 for the geometry of a global monopole with 
a = 0.5. 



6 Conclusion 

In the present paper we have considered the one-loop vacuum effects for a massive scalar field 
with general curvature coupling parameter on background of the (D + l)-dimensional global 
monopole with non-trivial core structure. The previous papers on the investigation of the 
vacuum polarization by the gravitational field of the global monopole are concerned with the 
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Figure 4: The renormalized energy density, a £>+1 (TQ ) ) ren inside the core for D = 3 massless 
scalar field as a function of f/aa in the flower-pot model with a = 0.5. The full/dashed curves 
correspond to minimally/conformally coupled scalars. 

idealized point-like model, where the curvature has singularity at the origin. The exception is 
the Ref. ^H|> where the vacuum densities for a massless scalar field are studied outside the 
monopole core with the interior de Sitter geometry. Here we consider the general spherically 
symmetric static model of the core with finite thickness, described by the line element P|). and 
investigate the vacuum properties in both exterior and interior regions. Among the most im- 
portant characteristics of these properties, which carry an information about the core structure, 
are the VEVs for the field square and the energy-momentum tensor. In order to obtain these 
expectation values we first construct the positive frequency Wightman function. In the region 
outside the core this function is presented as a sum of two distinct contributions. The first one 
corresponds to the Wightman function for the geometry of a point-like global monopole and 
the second one is induced by the non-trivial structure of the monopole's core. The latter is 
given by formula IJ29J) . where the tilted notation is defined by formula (|30|) with the coefficient 
from for the model without an infinitely thin spherical shell on the boundary of the core. 
This coefficient is determined by the radial part of the interior eigenfunctions and describes the 
influence of the core properties on the vacuum characteristics in the exterior region. In the case 
of the core model with a thin shell on the boundary the derivatives of the metric tensor compo- 
nents are discontinuous on the core surface. This leads to the delta function type contribution 
to the Ricci scalar and, hence to the equation for the radial eigenfunctions in the case of the 
non-minimally coupled scalar field. As a result, the radial eigenfunctions have a discontinuity 
in their slope at the core boundary. This leads to an additional term in the coefficient of the 
tilted notation which is proportional to the trace of the surface energy-momentum tensor (see 
Eq. ([SI])). 

By using the formula for the Wightman function, in section El we have investigated the 
influence of the non-trivial core structure on the VEVs of the field square and the energy- 
momentum tensor. As in the exterior region the local geometry is the same as that in the point- 
like global monopole model, the presence of the core does not lead to additional divergences 
for the points outside the core. As a result, the parts in these VEVs induced by the core are 
directly obtained from the corresponding part of the Wightman function for the case of the 
field square and by applying on this function a certain second-order differential operator and 
taking the coincidence limit for the energy-momentum tensor. These parts are given by formulae 
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(|33|) and (|39|) for the field square and the energy-momentum tensor respectively. They diverge 
as the boundary of the core is approached. The surface divergences in the VEVs of the local 
observables are well-known in quantum field theory with boundaries and are investigated for 
various boundary geometries. We have investigated the asymptotic behavior of the core induced 
VEVs near the core boundary and at large distances from the core. In particular, at large 
distances and for a massless scalar field with vq > 0, the ratio of the core induced and the 
point-like monopole parts decay as {r/a) 2v ° for the both field square and the energy-momentum 
tensor. For the special case with vq = this ratio decays logarithmically and long-range effects 
of the monopole's core appear. In the limit of strong gravitational fields corresponding to 
small values of the parameter cr, the behavior of the core induced parts is completely different 
for minimally and non-minimally coupled fields. The corresponding VEVs are suppressed by 
the factor exp [-(2/a)y/n(n + l)f ln(a/r)] for the non-minimally coupled scalar and behave like 
cr 1_D for the minimally coupled field. 

As an example of the application of the general results, in section we have considered 
a simple core model with a flat spacetime inside the core, so called flower-pot model. The 
corresponding surface energy-momentum tensor on the boundary of the core is obtained from 
the matching conditions and has the form given by Eq. Q56JI . The core induced parts of the 
exterior VEVs in this model are obtained from the general results by taking the function in 
the coefficient of the tilted notation form Eq. (|61|) . For the flower-pot model we have also 
investigated the vacuum densities inside the core. Though the spacetime geometry inside the 
core is Monkowskian, the non-trivial topology of the exterior region induces vacuum polarization 
effects in this region as well. In order to find the corresponding renormalized VEVs of the field 
square and the energy-momentum tensor we have derived a closed formula, Eq. (|72j). for the 
difference of the interior Wightman function and the Wightman function for the Minkowski 
spacetime. The subtracted function is finite in the coincidence limit and can be directly used 
for the evaluation of the VEVs of the field square and the energy-momentum tensor. The latter 
quantities are given by formulae l)74|) and (|79|) . As in the case of the exterior region, we have 
considered various limiting cases when the general formulae are simplified. In particular, we have 
shown that in the limit <r<l under the fixed value aa, which is the core radius for an internal 
Minkowskian observer, the renormalized VEVs of the field square and the energy-momentum 
tensor tend to finite limiting values. 
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A Contribution of bound states 

In this appendix we consider the contribution of possible bound states into the VEVs. For this 
states the quantity A is purely imaginary, A = irj, and the corresponding radial eigenfunction in 
the region r > a is the function A b ir~ n l 2 K Vl (j\r) with the normalization coefficient A^. To have 
a stable ground state we will assume that ry < m. From the continuity of the eigenfunctions at 
r = a one has 

Ri(a,i V )=A bl a- n / 2 K Ul ( V a), (82) 
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and from the continuity of the radial derivative we see that the possible bound states are solutions 
of the equation 

K ul ( V a)=0, (83) 
with the notation from (|30|L The normalization condition for the bound states is as follows: 



pa poo 

/ dr e- u+v+( - D -^ w Rf(r, i V ) + A bl a°- 1 / dr rK, 

Jro J a 



Sw-sskj- (84) 



from which the normalization constant can be found. In order to evaluate the integrals in this 
formula we note that for the solution f u i(r) to radial equation (jllj) the following formula takes 
place 



e -u+v+(D-l)w 



dre^+P-^faWrtir) = — 2 — [faW^r) - Li(r)fU r )} ■ ( 85 ) 



u>i — to 

In particular, in the limit uj\ — > u one finds 

dre - u+vH D-i) Wf 2 i{r) = e _ 



d d 
f u i(r)j^f u i(r) - f u i(r)—f ul (r) 



Applying to the integrals in Eq. (|84|1 this formula and using the continuity of the radial eigen- 
functions at r = a, for the normalization coefficient one finds 

A% = 7? " 1 " g/ -' M . (87) 

uoN{m k )(d/ drj)K Vl (rja) 

To obtain this formula we have used the relation 

K n ( V a) = l/I n {rja), (88) 

valid for the solutions of Eq. (|83|). As a result, for the contribution of the bound state with 
A = irj to the Wightman function we have the formula 



0.1— .D 00 21 4- n 
(v(xMx')) hs = -^^ 7 —^C? /2 (co S 8) 

1=0 

~& 

-K Vl {w)K Ul {w')- (89) 



nS D ^ (rr') n / 2 

X = 

(d/drj)K Ul (77a) ^Jm 2 - if 

In the case when several bound states are present the sum of their separate contributions should 
be taken. Now the Wightman function is the sum of the part coming from the modes with real 
A given by Eq. (|28[) and of the part coming from the bound states given by Eq. (|89|) . In order 
to transform the first part we again rotate the integration contour in Eq. (|28|) by the angle tt/2 
for s = 1 and by the angle —tt/2 for s = 2. But now we should take into account that the 
integrand has poles at A = ±zry which are zeroes of the functions H^' (Aa) in accordance with 
Eq. (|83|). Rotating the integration contour we will assume that the pole (— l) s irj, s = 1,2, on 

(s) 

the imaginary axis is avoided by the semicircle C p in the right half plane with small radius p 
and with the center at this pole. The integration over these semicircles will give an additional 
contribution 

_°Y^+JL C ^ 2 (cose)Y dXX 6 - 44^L^)(Ar)g(f)(ArQ. (90) 
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By evaluating the integrals in this formula it can be seen that this term exactly cancels the 
contribution (|89|) coming from the corresponding bound state (for the similar cancellation in 
the Casimir effect with Robin boundary condition see Ref. [27]). Hence, we conclude that the 
formulae given above for the core induced parts in the VEVs are valid in the case of the presence 
of bound states as well. 

In order to see the possibility for the appearance of bound states in the flower-pot model, we 
note that introducing a function Fi(r) = r^ D fi(r) with r* = r outside the core and r* = r 
inside the core, equation (jllj) for the radial part of the eigenfunctions is written in the form of 
the Schrodinger equation. The corresponding effective potential is equal [vf + (n — l)/4]/cr 2 r 2 
in the exterior region and [1(1 + n) + (n 2 — l)/4]/f 2 in the interior region. Under the conditions 
uf ^ and n > assumed earlier, the potential is non-negative and, hence, in the flower-pot 
model no bound states exist. 
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